PRESTRESS TECHNOLOGY 





LE EXER TT ee 


CH NZE. 》2015 年 第 5 期 总 第 112 期 





基于 群 论 的 预 应 力 结构 广义 特征 值 分 析 
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摘 要 :由 于 常规 自 振 分 析 方 法 未 能 充分 利用 结构 的 固有 对 称 性 ， 当 结构 自由 度 增加 时 ,计算 消耗 显著 提 
升 。 本 文 基于 群 论 方法 ， 提 出 一 种 分 析 对 称 型 预 应 力 结构 动力 特性 的 高 效 方法 。 首 先 ， 结 合 一 致 质量 矩阵 
及 切线 刚度 矩阵 建立 了 预 应 力 结构 的 广义 特征 方程 ， 并 考虑 初始 预 应 力 对 结构 的 影响 ， 以 求解 结构 的 自 振 
频率 及 振 型 。 随 后 ， 建 立 对 称 坐 标 系 ， 将 刚度 矩阵 及 质量 矩阵 分 解 为 一 系列 分 块 对 角 化 矩阵 。 由 于 各 分 块 
TERREI, 广义 特征 值 问题 的 求解 难度 显著 降低 ， 从 而 能 高 效 求解 结构 的 自 振 频 率 及 相应 振 型 。 数 
值 算 例 阐明 了 所 提出 方法 的 基本 计算 过 程 及 显著 优势 。 与 有 限 元 结果 及 常规 方法 所 得 结果 比较 后 可 知 ， 基 


于 群 论 的 对 称 方法 准确 、 高 效 。 
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现代 预 应 力 结构 主要 由 索 、 杆 、 梁 等 构件 组 
成 ， 在 世界 各 地 逐渐 得 以 应 用 ， 它 们 因 其 独特 的 
结构 外 形 以 及 较 强 的 跨越 能 力 受到 建筑 师 和 结构 
师 的 青睐 。 通 常 ， 在 初始 态 下 这 些 预 应 力 结构 包 
含 内 部 机 构 位 移 模 态 ， 当 施加 合适 的 预 应 力 后 ， 形 
成 一 定 的 结构 刚度 ， 从 而 达到 稳定 的 平衡 状态 。 

最 近 ， 预 应 力 结构 受 广泛 研究 ， 主 要 包括 形 
态 分 析 ( Pellegrino and Calladine 1986; Pellegrino 
1990) 、 稳 定性 分 析 (Gunnar 1999; Zhang et al. 
2009b ) 以 及 优化 分 析 ( Yuan and Dong 2003; 
Kitipornchai et al. 2005 ) 等 。 然 而 ， 针 对 现代 预 
应 力 结构 的 动力 分 析 研 究 却 极 少 。 实 际 上 ， 对 于 
一 个 结构 而 言 ， 固 有 振动 是 非常 重要 一 项 内 在 力 
学 特性 ， 它 直接 影响 结构 在 动态 作用 下 的 反应 。 
因此 ， 有 关 预 应 力 结构 固有 振动 特性 的 分 析 十 分 
重要 。 有 学 者 采用 有 限 元 方法 分 析 了 Giger HRA 
顶 的 自 振 特性 ( Luo and Wang 2005 ) ， 研 究 已 经 
发 现 该 类 型 结构 的 自 振 频率 相当 低 ， 并 且 分 布 集 
"P. Zhang 等 人 (Zhang 2007 ) 对 张 弦 梁 结构 的 固 
有 频率 进行 了 研究 ， 文 中 解释 说 明了 索 结构 的 张 
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拉 以 及 不 同 边界 条 件 对 动力 特性 的 显著 影 啊 。 
Wu 和 Sasaki ( Wu and Sasaki 2007 ) 采用 数值 方法 
对 索 拱 结构 的 静态 和 动态 特性 进行 了 研究 ， 他 
们 所 求 得 的 拱 的 固有 振动 的 结论 与 试验 值 吻合 
非常 好 。 

”从 数学 角度 来 看 ， 结 构 的 目 振 频 率 可 以 通过 
求解 关于 结构 切线 刚度 矩阵 和 质量 矩阵 的 广义 特 
征 值 方程 直接 获得 。 然 而 ， 由 于 结构 目 由 度 的 增 
加 ， 计 算 需 求 也 显著 增加 。 在 这 种 情况 下， 我 们 
必须 采取 有 效 的 求解 方法 。Kaveh 和 Rahami 
(2007a, b ) 研究 了 三 对 角 线 矩阵 和 五 对 角 线 和 矩 
阵 ， 并 且 提 出 了 有 效 的 分 解 方法 。 之 后 ， 有 学 者 
提出 了 结合 拓扑 学 和 图 论 的 数值 解法 ， 用 于 解决 
圆柱 形 网 格 结构 的 特征 值 问 题 ( Kaveh and Rahami 
2008) 。 由 于 大 部 分 结构 要 么 是 对 称 的 ， 要么 是 
由 重复 的 子 结构 组 成 ， 群 论 可 有 效 简化 这 些 对 称 
结构 的 广义 特征 值 分 析 。 在 物理 和 化 学 领域 ， 基 
于 群 论 的 对 称 学 方法 已 经 被 证 明 是 一 个 非常 强大 
的 数学 工具 ( Ceulemans and Fowler 1991 ) , 然 
而 ， 在 土木 工程 领域 却 鲜 有 应 用 。 最 近 几 年 ， 群 
论 逐 渐 被 用 于 一 些 结构 力学 的 研究 ( Healey and 
Treacy 1991; Zingoni 1996, 2005, 2008; Mohan 
and Pratap 2004; Kaveh and Nikbakht 2007, 2010; 
Kaveh and Rahami 2010a ) o Kangwai FA ( Kangai 
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et al. 1999 ) 具体 解释 了 如 何 运 用 群 论 求解 对 称 
结构 的 对 称 坐 标 系 。 最 近 ， 也 有 文 著 综 述 了 群 
论 方法 在 结构 工程 分 析 研 究 中 的 应 用 情况 和 优 
Ñ, ( Zlokovic 1989; Zingoni 2002, 2009) 。 


然而 ， 很 少 有 研究 涉及 到 对 称 型 预 应 力 结 


构 的 求解 ， 这 是 由 于 预 应 力 结构 借助 合理 的 预 
应 力 来 获得 其 自身 稳定 。Pandia Raj 和 Guest 
( Pandia Raj and Guest 2006 ) 利用 群 论 来 简化 张 
拉 整 体 结构 的 找 形 过 程 。 他 们 发 现 ， 平 衡 矩 阵 
和 应 力矩 阵 均 能 够 被 对 角 化 。Kaveh 等 人 《Kaveh 
et al. 2010 ) 利用 图 论 进 了 预 应 力 索 网 结构 的 频 
率 分 析 ， 其 中 相关 和 矩阵 可 实现 完全 的 分 块 对 角 
化 ， 因 而 能 高 效 地 获得 自 振 频 率 。 在 此 研究 基 
础 上 ， 本 文 将 对 具有 复杂 对 称 性 的 现代 预 应 力 
结构 进行 更 为 系统 的 广义 特征 值 分 析 。 
2 基本 理论 
2. 1 基本 假定 及 广义 特征 值 方程 
需 指出 ， 文 中 采取 以 下 假定 : 
(a) BUR Tr SUA SB HE 
(b) RAE HHMI, HEH- 
应 变 关 系 为 线 弹 性 ; 
(c) 各 单元 自重 平均 分 配 到 两 端 ; 
Cd) 结构 的 固有 振动 模式 为 简 谐 振动 。 

忽略 阻尼 作用 ， 预 应 力 结构 的 广义 特征 值 
方程 为 : | 

K$-w Mó-F-0 (1) 
式 中 ，MM 为 具有 m 个 节点 的 d 维 结构 (d=2 或 
d=3 ) 的 md x md 对 称 质量 矩阵 ; 

有 为 md x md 对 称 切线 刚度 矩阵 ; 

w 为 结构 自 振 频率 构成 的 向 量 ; 

fb 为 相应 的 Pad x 1 振 型 矩阵 。 

具有 2 个 单元 的 预 应 力 结构 的 切线 刚度 矩阵 
可 以 与 成 : 

K=HGH' + (CT tC) QL (2) 
式 中 ， 五 为 md x bY Se, HRCA HS 
点 外 力 的 关系 ; 

G 为 b x 5 对 角 和 矩阵 ， 其 包含 单元 的 修正 轴 向 
刚度 ( Guest 2006; Zhang et al. 2009a ) ; 

IL 为 d xd 的 单位 矩阵 ; 
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GO 表示 克 罗 内 克 积 ; 

t= t/ [是 b x 5 的 应 力矩 阵 ， 它 反应 了 单位 长 
度 的 应 力 等 级 !， CAKE, 

C 是 b x m 的 拓扑 矩阵 ， 它 描述 了 几何 拓扑 天 
系 ， 如 果 单 元 连接 节点 i 和 j ， 那么 矩阵 的 第 
i 行 、j 行 、k 行 分 别 为 1、-1 和 0。 

2.2 群 论 及 其 矩阵 表示 

采用 群 论 方法 (Kaveh and Nikbakht 2007; 
Zingoni 2009 ) 或 图 论 - 代 数 方法 (Kaveh and Fazli 
2010; Kaveh and Raham 2010b ) ， 能 够 对 对 称 结 
构 进 行 深入 系统 的 研究 。 本 市， 首先 将 对 群 论 及 
其 矩阵 表示 的 概念 作 一 些 简单 介绍 。 

群集 C=[fg, i=1,2,..., n) 是 系列 独立 元 素 的 集 
合 ， 并 满足 以 下 四 项 准则 | 

(1) 必须 包含 单位 元 素 : 存在 ECG， 对 于 
任何 元 素 g CG, Hef Eg - E-E + gg; 

(i) 任何 两 个 元 素 的 运算 结果 仍 是 群集 的 
TH, Pile, CG, gj CG, 那么 gj - 2, CG; 

(ii) 每 个 元 素 在 群集 中 都 有 其 逆 元 素 ， 如 
gjC G， 则 g-1C G， 反 之 亦 然 ; | 

(iv) 群集 元 素 遵 循 结合 律 ， 例 如 : (ge 
gz ) ”83=81 7 (g ] g3 ) o 

为 了 描述 一 个 结构 的 对 称 性 ， 我 们 将 各 个 对 
称 操作 视 为 群集 的 独立 元 素 。 对 称 操作 是 一 种 特 
定 的 操作 ， 在 对 称 操作 下 ， 结 构 变 换 后 的 几何 构 
型 与 原 几何 构 型 完全 等 效 ( 几何 上 无 法 识别 ) 。 
我 们 可 以 用 和 矩阵 的 方式 描述 在 对 称 操作 下 节点 坐 
标 X 的 变换 和 节点 荷载 F 的 变换 : 

X-RX, HF-RF (3) 
式 中 , 天 和 天 分 别 是 在 对 称 坐 标 系 下 节点 坐标 矩 
阵 和 节点 和 荷载 矩阵 ， 和 矩阵 R 是 对 应 的 转换 矩阵 。 
在 对 称 操作 作用 下 ， 转 换 和 矩阵 R 可 表示 为 : 

R=S, GOr (4) 
式 中 ，5, 是 描述 一 个 结构 固有 对 称 性 m x m 的 拓 
扑 和 矩阵。 在 对 称 操作 下 ， 乔 将 笛 卡 尔 坐 标 系 中 某 
qaa RET ui, 则 拓扑 矩阵 可 通过 取 第 /7 íT 
E58 i 行 分 别 为 1 和 0 获得 。7 是 从 节点 7 坐标 转换 
到 节点 i 坐标 的 d x 4 的 局 部 转换 矩阵 。 在 此 ， 以 
一 个 简单 的 二 维 对 称 索 杆 结构 为 例 ， 如 图 1 所 








示 。 在 绕 原点 0 旋转 x 角度 的 对 称 操作 下 ( 用 C，- 


表示 ) ， 节 点 1 到 节点 6 分 别 由 初始 节点 3、4、 
1、2、6 和 5 转换 而 来 。 因 此 ， 对 这 一 特定 的 对 称 
操作 ， 转 换 矩 阵 表 示 为 : 


001000 (5) 
000100 

peli 9 9 9 0 aoe 
0 10 00 0 sin(1) cos( 7 ) 
000001 
00001 0 





图 1 二 维 C>, 对 称 索 杆 结构 


在 对 称 操作 (C,) F, AR (5) 中 2x2 的 
局 部 转换 矩阵 使 得 每 个 节点 的 x RI y 轴 变 成 相 
反 的 方向 。 实 际 上 ， 所 有 类 似 于 公式 C5) 中 的 
转换 矩阵 均 为 稀世 和 矩阵 ， 因 此 能 够 被 继续 化 简 , 
其 分 解 后 的 子 矩 阵 被 称 为 不 可 约 表示 Tr“。 当 
然 ， 对 于 不 同 的 点 对 称 群 ， 所 对 应 的 对 称 操作 、 
不 可 约 表示 均 能 够 从 一 些 群 论 书籍 中 直接 查阅 得 
到 ， 例 如 人 参考 文献 Altmann and Herzig (1994 ) 和 
Kettle ( 1995) 。 

在 土木 工程 领域 ， 传 统 的 对 称 预 应 力 结构 包 
括 周期 对 称 性 结构 ( 如 重复 的 张 拉 整体 结构 ) , 
或 环 向 轴 对 称 〈 如 索 穹顶 结构 ) 。 需 指出 ， 周 期 
对 称 结构 中 包含 重复 的 基本 单元 ， 在 同一 方向 以 
同等 距离 平移 对 称 操作 ， 移 动 结构 所 有 的 节点 和 
单元 。 为 了 使 用 平移 对 称 性 研究 周期 对 称 结构 ， 
需要 忽略 影响 局 部 力学 性 能 的 结构 边界 条 件 。 值 
得 一 提 的 是 ， 本 文 仅 考 虑 属于 环 辐 对 称 群 C 的 
对 称 预 应 力 结构 。 三 种 不 同 的 对 称 操作 构建 了 环 
[HD ERHEC, : 全 等 操作 EE， 沿 主轴 逆 时 针 旋 转 操 
WC? (i=1,2…n-1) ， 以 及 沿 包 含 相 应 主轴 的 
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平面 镜像 操作 o，(j=1, 2…n ) 。 这 2n 个 对 称 操作 
满足 群集 的 所 有 四 个 准则 要 求 ， 因 此 能 够 构建 成 
一 个 封闭 的 对 称 群集 。 文 中 ， 以 笛 卡 尔 坐 标 系 中 
的 z 轴 规定 为 主轴 ， 而 平面 xoy 与 用 于 镜像 操作 的 
镜面 平面 相互 垂直 。 | 

表 1 归 纳 了 环 向 对 称 群 C,, 的 不 可 约 矩 阵 表 
示 。 对 称 群 C 总 计 有 2a+2 个 一 维 的 不 可 约 矩 阵 
表示 及 p 个 二 维 的 不 可 约 矩 阵 表示 ， 其 中 


n-1-o _{0 when7 is odd 
= 7 HOF], when M is even (6) 


Ri 0 环 向 对 称 群 C,, 的 不 可 约 和 矩阵 表示 工 


pub .w M poo ea 位 移 5 转角 

A, 1 1 z 

A; 1 1 «i R, 

B, 1 iy (-1)/ 

B, 1 Cay" (-1)” 

E, i -G y G uS (x, y) 
01 S G S -C  (R.R) 

E, i ues. Cie. Si 

-2,"s,p 0 1 Six Cx Six -Cp 


ik: PB E BSBÁTAA BAUR PRIN ow, Ae 
Sj, sin e) » Cy=cos (um 


2.3 对 称 型 刚度 矩阵 及 对 称 型 质量 矩阵 

” ”通过 对 群 论 中 非常 重要 的 广义 正 交 定理 ， 可 
以 推导 出 关于 一 维 不 可 约 和 矩阵 表示 T “的 对 称 子 
空间 ,为 : 

Vik CX ry +R) (7) 
sth p“sA,, A, B,, By, XTPRPRVEISE, C, o 
C, , Ois 05, Oms JERE PRB (x ) 用 于 求 
ILE et FB BE x 的 列 空间 。 

此 外 ， 关 于 二 维 不 可 约 和 矩阵 表示 的 对 称 子 空 
间 为 : 

Vien sFOE TIG, R) (8) 
式 中 j=1，2 ， 且 对 称 操作 与 式 (7 ) 中 相同 。 因 
此 ， 用 于 将 坐标 和 荷载 转换 到 对 称 坐 标 系 下 的 正 
AEHHRORWEVON : 

ym |^ CV. (9) 

随后 ， 可 将 对 称 型 坐标 矩阵 和 外 力矩 阵 写 成 


”如 下 形式 : . 
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X-V'X, F=V'F (10) 
考虑 到 位 移 矩 阵 6 可 表示 为 节点 坐标 的 线性 
变换 ， 可 以 得 出 : 


b=V 0 (11) 
因此 ， 求 得 对 称 型 刚度 矩阵 为 : 
Ky es 
K, 0 
K-V'KV- Ke, (12) 
0 K, 


根据 本 文 第 2.1 节 所 列 的 基本 假定 ， 可 知 单元 
和 节点 的 目 重 己 对 应 的 节点 相关 。 所 以 ， 与 节操 
位 移 的 转换 矩阵 类 似 ， 可 将 质量 矩阵 转换 成 对 称 
型 分 解 矩阵 ， 即 : 


( 13) 


式 C12) 和 式 C13) 表明 ， 对 称 型 刚度 矩阵 
有 和 对 称 型 质量 矩阵 好 已 经 被 分 解 为 分 块 对 角 化 
矩阵 ， 并 且 包 含 了 p+2a+2 个 沿 对 角 线 分 布 的 子 
矩阵 。 因 此， 可 通过 独立 求解 每 个 子 矩 阵 玉 及 其 
相应 的 M, ， 求 解 相 关 和 矩阵 的 特征 值 ， 即 : 
v, M, $.-K, 6, ( 14) 
式 中 i=4|，A，,，.….，E,。 因 而 ,我们 可 以 有 效 获 
得 属于 对 称 子 空间 V (i ) 的 自 振 频 率 o A He 
态 办 。 由 于 相似 变换 并 不 改变 矩阵 的 特征 值 ， 所 
以 求 得 的 特征 值 即 为 原 结构 的 自 振 频率 ， 也 就 是 
抽 =w。 原 结构 的 自 振 模 态 可 表述 为 : 

p=V (i) d, (15) 

以 上 分 析 过 程 表明 ， 原 有 的 关于 3m 个 自由 
度 的 结构 自 振 分 析 可 以 通过 独立 求解 p+2a+2 个 广 
义 特征 方程 获得 。 
3 数值 算 例 

为 了 验证 本 文 所 提 对 称 方法 的 准确 性 和 有 效 
性 ， 本 节 将 列 出 一 些 关于 对 称 预 应 力 结构 的 实 





CPRESTRESS TECHNOLOGY ) 
LE CFEC Ee ne 


(38.2 324 X. 》2015 年 第 5 期 总 第 112 期 





例 。 其 中 ， 结 构 构 型 〈( 主要 是 节点 和 单元 的 编 
号 ) 将 合理 排列 ， 从 而 满足 对 称 法 的 要 求 ， 并 且 
便于 生成 对 称 拓扑 短 阵 (公式 (4) 所 定义 的 
S.) 。 此 外 ,利用 传统 方法 进行 了 分 析 和 比较 ， 
从 而 验证 利用 群 论 方法 所 得 结果 的 可 靠 性 。 应 当 
指出 ， 对 称 法 和 传统 数值 法 均 运行 于 相同 的 硬件 
平台 (英特尔 奔腾 双核 2.0GHZ ) 。 在 此 还 需 指 
出 ， 算 例 中 除非 特殊 说 明 ， 索 和 撑 杆 的 弹性 模 量 
分 别 为 Ew。=1.9e5MPa , E,,,-2e5MPa ， 所 有 单 
元 的 质量 密度 均 为 p =7850kg/m’。 
3.1 C,, 对称 索 杆 结构 

首先 ， 以 一 个 简单 的 二 维 C,, 对 称 索 杆 结构 
为 例 ， 详 细 描 述 采用 群 论 方法 求解 对 称 预 应 力 结 
构 自 振 频 率 的 基本 过 程 。 结 构 如 图 1 所 示 ， 图 中 
给 出 了 节点 和 单元 的 编号 以 及 边界 条 件 。 水 平 索 
(单元 1-2 ) 和 撑 杆 (单元 7-8 ) 的 松弛 长 度 为 
Ly-L;-2000mm, SHR ( 单元 3-6 ) 的 松弛 长 度 
为 上 ,=1414.21mm。 索 的 截面 面积 ， 
4. =500mm2， 撑 杆 的 截面 面积 ; 
Å =3000mm2。 我 们 对 水 平 索 施加 初始 预 应 力 大 
小 为 100kN ， 其 它 单元 的 预 应 力 大 小 可 根据 来 自 
平衡 矩阵 H 的 自 应 力 模 态 来 确定 。 

由 节点 1-8 构 建成 的 C,, 对称 索 杆 结构 在 以 下 
对 称 操作 作用 下 能 转换 成 等 效 的 构 型 : 全 等 变换 
CE) 、 绕 原点 旋转 x (C) 、 对 x WAR Cou) 
以 及 对 》 轴 镜像 (o5) 。 四 个 独立 对 称 操作 { EE， 
C,，o1，0;] 组 成 了 对 称 群 C,,。 如 表 1 ( 取 n=2 ) 
和 公式 (9 ) Bros, TEC, 对称 群 中 的 不 可 约 表 
示 ，C,, 对 称 索 杆 结 构 的 正 交 转换 矩阵 V 可 推导 求 
f$, WX (16)。 





A A B, B, 

OO TG 3109. 

1 0:1 0'1 0-10 
V=0.5x 0 1:0 -110 -1:0 1 

a 0i- 00-1 0-1 0| ‘16? 

0 1:0 1:0 -1:0 = 

1 O:-1 0:1 0-1 0 

0 -110 1'0 -110 1 

10:1 0-1 0-1 0 





EEV 中 含有 4 个 对 称 子 空间 ， 其 中 在 
V CA) 中 ， 和 荷载 和 位 移 完 全 对 称 ; V CAD 中 ， 
它们 为 C, 旋转 对 称 ; V (B) 中 对 称 操作 oa, 和 
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V(B,) 中 o, 被 继续 保留 。 由 此 可 知 ， 对 称 型 刚 
度 和 矩阵 和 对 称 型 质量 矩阵 可 通过 公式 (2 ) 、 
(12-13) 和 (16 ) 计算 求 得 ， 其 子 和 矩阵 如 表 2 
所 列 。 


表 2 ” 子 答 阵 的 对 称 匹 配 刚 度 矩 阵 以 及 质量 矩阵 


FIERE A, A; B, B, 
E E, able Acable E hi p hi » v bi 07 | 
.(x-——————— 
if Li 07 27 ||-07 07 j1-07 0.7 || 07 27 


M, xe A uu (Lg +Ly) E | in d [s 0 | 区 d 
| 0 0.5] 10 05} [005][ 005 


分 块 对 角 化 刚度 矩阵 和 质量 矩阵 M 包含 
四 个 沿 对 角 线 独立 分 布 的 2 x 2-F BE. Piu, R 
们 可 以 通过 求解 子 矩 阵 K AUR, (i= A,, 
A,, B,, B,) 来 求解 特征 值 ， 具 体 公 式 如 下 : 

w M,ó-K,d;-05X i2 A,, A, B,, B, (17) 

应 当 指出 ， 原 结构 共 具 有 八 个 自由 度 ， 其 频 
率 分 析 被 四 个 独立 的 二 维 问 题 所 替代 。 表 3 详 
细 列 出 了 有 预 应 力 结构 和 无 预 应 力 结构 的 计算 
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由 表 3 可 知 ， 对 称 法 求解 结果 与 数值 法 完 
全 一 致 。 与 此 同时 ， 对 称 法 计算 所 耗 时 间 仅 为 
0.004s， 而 传统 数值 法 所 耗 时 间 高 达 0.01S6s 。 从 
数学 的 角度 来 讲 ， 群 论 方法 与 传统 数值 法 一 样 ， 


“但 由 于 相关 矩阵 已 被 分 块 对 角 化 ， 能 够 被 独立 求 


解 ， 因 而 群 论 方法 更 为 高 效 。 

对 模型 A 和 型 8 进行 比较 可 以 发 现 ， 由 于 第 1 
阶 自 振 频 率 完全 依赖 于 预 应 力 的 等 级 ， 所 以 C，, 
对 称 索 杆 结构 可 通过 预 应 力 实现 稳定 。 由 于 存在 
预 应 力 ， 结 构 的 低 阶 自 振 频 率 略 增加 了 一 些 。 

图 2 描绘 了 结构 所 有 的 自 振 振 型 ， 并 根据 式 

(15) 给 出 了 各 振 型 所 属 的 对 称 子 空间 。 图 中 短 

箭头 表示 各 节点 的 相对 运动 趋势 ， 可 以 发 现 结构 
振 型 与 其 对 称 子 空间 关联 的 对 称 特 性 一 致 。 模 态 
5 和 8 具有 全 对 称 性 ， 这 是 由 于 他 们 来 自 对 称 子 空 
间 4,。 类 似 地 ， 模 态 1 和 4 具有 旋转 对 称 性 ， 模 态 
3 和 7 关于 x 轴 对 称 ， 而 模 态 2 和 6 关于 y 轴 对 称 。 


表 3 对 称 索 杆 结构 的 自 振 频 率 


模 态 1 2 3 
对 称 法 9.14 138.72 163.10 
模型 A 称 
数值 法 9.14 138.72 163.10 
对 称 法 0 138.51 162.95 
模型 B TERIA 
数值 法 0 138.51 162.95 


TE: 模型 A 与 模型 B 分 别 定义 为 有 预 应 力 与 无 预 应 力 结构 。 





Mode 5, subspace A, Mode 6, subspace B, 


4 5 6 7 8 
237.17 325.36 341.40 729.54 729.68 
237.17 325.36 341.40 729.54 729.68 
237.17 325.29 341.48 729.46 729.66 
237.17 325.29 341.48 729.46 729.66 


Mode 4, subspace B, 





Mode 7, subspace B, Mode 8, subspace A, 


图 2 对 称 索 杆 结构 的 自 振 振 开 


— @ -—— 
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3.2 Co 对 称 索 网 结构 

本 节 将 分 析 一 个 由 12 个 节点 和 12 个 索 单 元 组 
成 的 三 维 索 网 结构 ， 结 构 具 有 Ce, 对 称 性 。 图 3 给 
出 了 结构 的 节点 和 单元 编号 ， 并 给 出 了 由 6 个 主轴 
线 1-6 (o-o) 定义 的 对 称 边界 条 件 。 在 此 ， 取 
径 向 索 (单元 1-6 ) 的 松弛 长 度 为 Lr=4000mm , 
Rie) (单元 7-12 ) 的 松弛 长 度 为 L,=3000mm 。 
索 的 横 截 面 面 积 均 为 4.upe=1963.44mm2， 所 有 单 
元 的 初始 预 应 力 大 小 均 为 566kN 。 


NG5 y 
\ 





图 3 三 维 C6v 对 称 索 网 结构 
对 称 群 Co 由 12 个 独立 的 对 称 操作 组 成 ， 即 
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IRELE, Cl, C2, CÈ, Co, Co, On Os 04, 
O, Os, Og])o WIRI (u6) MAÑ (9) 的 
不 可 约 表示 ，Ca 对 称 索 网 结构 的 正 交 转换 矩阵 
VA: 

V=(V(A,), WA), WB)), VB), WE) 
V( Ej), VE), WEn)] (18) 

其 中 ， 各 不 可 约 表示 关联 的 对 称 子 空间 如 
表 4 所 示 。 

由 表 4 可 知 ， 正 交 转 换 矩 阵 Y 中 存在 8 个 对 称 
子 空间 。 节 点 荷载 和 位 移 在 V (4) ) 中 是 完全 对 
称 的 ， 并 在 V (A,) 中 是 旋转 对 称 的。 在 V (Bj) 
中 沿 轴 1、3、5 保持 镜像 对 称 及 C, 旋 转 对 称 ; 在 
V CB,) 中 沿 轴 2、4、6 镜像 对 称 及 C, 旋 转 对 称 。 
uo, SEV CE, 中 保持 o, 对 称 ， 在 
V (En) 中 保持 o ,对称 ; 它们 在 子 空间 V CES) 
中 沿 轴 线 1、4 保 持 C, 旋 转 对 称 ， 在 V(E,, ) 中 亦 
保持 C, 旋 转 对 称 。 | 

由 于 所 有 非 约 束 节点 具有 旋转 对 称 属性 ， 结 
构 的 对 称 型 质量 矩阵 仍然 保持 原 对 角质 量 和 矩阵 不 
变 ， 且 矩阵 中 元 素 M -0.5pA uu, CLg2L,) 。 表 5 
给 出 了 Cs, Mt PR PI Ba PY MT HR DE AE PA EE E E 
子 块 : 


X4 — C， 对 称 索 网 结构 的 不 同 对 称 子 空间 


对 称 子 空间 矩阵 表示 

WA) [-0.41, 0, 0, -0.20, 0.35, 0, 0.20, 0.35,0, 0.41, 0, 0, 0.20, -0.35, 0, -020, -035, OY 

VCA;) [0, -0.41, 0, ~0.35, -0.20, 0, ~0.35, 0.20, 0, 0, 0.41, 0, 0.35, 0.20, 0, 0.35, -0.20, 0; 
0, 0, -0.41, 0, 0, -0.41, 0, 0, -0.41, 0, 0, -0.41, 0, 0, -041, 0, 0, -0.41]" 

VB; ) [-041, 0, 0, 0.20, -0.35, 0, 0.20, 0.35, 0, -0.41, 0, 0, 0.20, -035, 0, 0.20, 035, 0J 

V(B,) [0, -0.41, 0, 0.35, 0.20, 0, -0.35, 0.20, 0, 0, -0.41, 0, 0.35, 020, 0, -0.35, 0.20, 0; 
0, 0, -0.41, 0, 0, 0.41, 0, 0, -0.41, 0, 0, 0.41, 0, 0, -0.41, 0, 0, 041 f 

VE) [0, 0, 0, -0.43, -0.25, 0, -043, 025, 0, 0, 0, 0, -0.43, -025, 0, -0.43, 0.25, 0; 


0, 0, 0, 0, 0, 0.50, 0, 0, 0.50, 0, 0, 0, 0, 0, -0.50, 0, 0, -0.50; 

-0.58, 0, 0, -0.14, 0.25, 0, -0.14, -0.25, 0, -0.58, 0, 0, -0.14, 0.25, 0, -0.14, -025, 0]^ 
WE) [0. 0, 0, 0.25, -0.43, 0, -025, -0.43, 0, 0, 0, 0, 0.25, -043, 0, -025, .043, 0; 

0, -0.58, 0, -0.25, -0.14, 0, 0.25, -0.14, 0, 0, -0.58, 0, -0.25, -0.14, 0, 025, -0.14, 0; 

0, 0, -0.58, 0, 0, -029, 0, 0, 029, 0, 0, 0.58, 0, 0, 029, 0, 0, -0.29 J 
VC Ex, ) [0, 0, 0, 0.43, 0.25, 0, -0.43, 025, 0, 0, 0, 0, -0.43, -025, 0, 0.43, -025, 0; 

0, 0, 0, 0, 0, 0.50, 0, 0, -0.50, 0, 0, 0, 0, 0, 0.50, 0, 0, 0.50; 

0.58, 0, 0, -0.14, 0.25, 0, 0.14, 0.25, 0, -0.58, 0, 0, 0.14, -0.25, 0, -0.14, -0.25, OF 
Vi Ej) [0, 0, 0, -0.25, 0.43, 0, -0.25, -0.43, 0, 0, 0, 0, 025, ~0.43, 0, 0.25, 0.43, 0; 

0, 0.58, 0, -025, -0.14, 0, -0.25, 0.14, 0, 0, -0.58, 0, 025, 0.14, 0, 025, -0.14, 0; 

0, 0, -0.58, 0, 0, 029, 0, 0, 029, 0, 0, -0.58, 0, 0, 029, 0, 0, 029]" 


一 一 和 -一 一 








表 5 
FFE A, A> 
ur E do Á cable | -0.03 0 | 
K, (x 
OC oL. >? [23.33] 0 -003 
TERE E, 

20 -10 0 

EN Eye A 

K ( x cable “cable ) -10 9.97 0 

100L, 

0 0 —0.03 





由 表 5 可 知 ， 对 称 型 刚度 矩阵 K 被 分 块 对 角 
化 ， 沿 着 主 对 角 方向 包含 两 个 1 x 1 的 子 和 矩阵 块 
Kan Kp; 两 个 2 x 2 的 子 矩 阵 块 局、Kp; 四 个 3 
x 3 的 子 和 矩阵 块 K; ~ Ken Ke, 和 Ks 。 所 以 特征 
值 可 以 通过 求解 每 个 子 矩 阵 K, 及 其 相应 守 ，( i= 
ur A), B,, B,, Ej, Ey, E, E) Ei 
作为 相互 独立 的 广义 特征 值 问题 来 求解 。 随 后 ， 
我 们 可 以 得 到 对 称 子 空间 V (i ) 的 频率 o, C ELS 
BU p, 。 自 振 模 态 可 表示 为 : 

9=V (i) 办， 其 中 i= (i=A,, A, B, B,, 
Ej. Ey, Ey, E) (19) 

图 4 给 出 了 所 得 结果 与 常规 数值 结果 的 对 比 
情况 。 可 知 ， 由 群 论 方法 所 求 得 的 频率 与 常规 法 
十 分 吻合 ， 二 者 之 间 的 最 大 误差 为 9.4e-8， 并 且 
出 现在 自 振 模 态 的 第 14 阶 和 18 阶 。 需 指出 ， 两 种 
方法 在 保持 同等 精度 的 情况 下 ， 常 规 数值 法 所 需 
计算 时 间 比 对 称 法 六 倍 还 多 ， 高 达 0.0781s 。 

34 


^  symmetry,with prestress 


30 numerical method, with prestress 
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图 4 对 称 索 网 结构 的 低 阶 自 振 频率 
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对 称 型 刚度 和 矩 阵 的 分 块 子 和 矩阵 


B, B, Ey 
Ho] | 39.97 0 | 9.97 0 -10 
0 -0.03 0 
0 —0.03 d 
-10 0 20 
E; E» 
29.97 Ü -10 13.33 ~10 0 
0 -003 0 -10 2997 0 
~10 O0 -13.33 0 0 —0.03 











由 于 该 索 网 结构 的 整体 平衡 矩阵 H 存 在 奇 晃 
值 ， 前 七 阶 自 振 模 态 对 应 着 七 个 一 阶 无 穷 小 机 构 
人 位移。 因此， 如 公式 (2) 所 列 的 切线 刚度 矩阵 
所 体现 的 ， 前 七 阶 自 振 频 率 对 结构 初始 预 应 力 较 
为 敏感 。 当 结构 不 存在 初始 预 应 力 时 ， 相 应 的 自 
振 频 率 为 零 ， 如 图 4 中 下 部 虚线 所 示 。 即 便 如 
此 ， 采 用 对 称 法 仍然 准确 求 得 了 无 预 应 力 索 网 
结构 的 所 有 结果 ， 与 常规 数值 法 的 最 大 误差 为 
2.57e -7。 

图 5 为 该 对 称 索 网 结构 的 前 12 阶 自 振 模 态 ， 
可 以 看 出 ， 结 构 的 第 1 阶 自 振 模 态 为 Ce, 对 称 ， 第 
3 阶 模 态 为 C¢ 旋转 对 称 ， 第 7 阶 和 第 12 阶 模 态 均 为 
C, 对 称 。 因 此 可 以 确定 ， 图 5 所 示 的 结构 自 振 模 
态 均 满足 相应 对 称 子 空间 的 特征 属性 。 

由 于 与 二 维 不 可 约 表示 关联 的 对 称 子 空间 总 
存在 重 根 ， 我 们 在 数值 计算 之 前 即 可 据 此 推测 出 
对 称 子 空 间 E,, 和 E,，(i =1，2 ) 将 会 产生 相同 的 
自 振 频 率 。 事 实 也 确实 如 此 ， 图 4 和 图 5 所 示 的 结 
果 证 明了 以 上 推断 。 模 态 2 和 模 态 4 具有 相同 的 自 
振 频 率 ， 模 态 8 和 模 态 9 具有 相同 的 自 振 频 率 ， 这 
是 由 于 它们 均 来 自 对 称 子 空间 El; 类 似 地 ， 模 态 
5 和 模 态 6 具有 相同 自 振 频率 ， 模 态 10 和 模 态 11 具 
有 相同 自 振 频率 ， 它 们 都 来 自 于 对 称 子 空间 Ek,。 
3.3 C, 对 称 Levy 型 索 穹顶 结构 

本 节 将 围绕 一 个 几何 构 型 更 为 复杂 的 结构 一 
BIC a 对 称 Levy 型 索 穹 项 结构 ， 进 行 其 广义 特征 
值 分 析 。Levy 型 索 帘 项 是 一 种 新 型 的 预 应 力 索 杆 
结构 ， 由 美国 工程 师 Levy 提 出 ( Gerardo 和 Levy 
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1992) 。 此 后 ， 在 世界 范围 内 (尤其 是 美 、 日 、 
韩 等 发 达 国 家 ) 已 相继 建成 了 一 些 典 型 的 Levy 型 
索 穹 顶 结 构 ， 例 如 美国 亚特兰大 的 乔治 亚 穹 项 ， 
结构 跨度 为 233.5m x 186m. 

图 6 给 出 了 一 个 Levy 型 索 容 顶 结 构 的 几何 构 
型 ， 该 结构 的 加 直径 为 100m ， 图 中 也 给 出 了 结 
构 的 几何 边界 条 件 和 对 称 性 。 该 Ci 对 称 Levy 型 
索 穹顶 结构 由 84 个 贸 接 节点 和 204 个 单元 组 成 ， 





Mode 1, subspace A, Mode 2, subspace Ej; 


Mode 10, subspace Ey 
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其 中 12 个 外 边界 节点 沿 三 个 方向 受到 位 移 约束 。 
所 有 单元 包括 3 种 撑 杆 和 10 种 索 ， 每 一 类 型 单元 
都 包含 12 个 绕 z 轴 环 加 分 布 的 单元 。 就 群 论 所 描 
述 的 对 称 性 而 言 ， 该 结构 拥有 1 个 全 等 变换 操作 
E，11 个 旋转 对 称 操作 Ci,- Ci,， 以 及 12 个 镜像 对 
称 操作 oj-o,， 如 图 6 (a) 所 示 。 

表 6 具 体 给 出 了 不 同类 型 单元 的 横 截 面 面 积 
和 初始 预 应 力 大 小 。 





Mode 11, subspace Æ» 


图 5 对称 索 网 结构 的 前 12 阶 月 振 振 型 





全 6 C15, Xt PK Levy WRS Til: Ca) 边界 条 件 及 对 称 条 件 ; 


(b) 三 维 视图 及 节点 编号 ; 





(c) 镜面 o-oo 中 的 不 同类 型 与 单元 长 度 ( mm ) 
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表 6 ”不 同类 型 单元 的 截面 积 与 初始 预 应 力 


单元 类 型 jl J2 J3 X1 X2 X3 
面积 /cm 36.7 32.7 24.9 42 36.7 16.7 


预 应 力 /kN 1343 644 392 901 688 248 


取 表 1 中 n=12 ， 并 结合 公式 (7-9 ) ， 可 以 
求 得 该 C,， 对 称 索 容 顶 结构 的 正 交 转换 矩阵 V。 
矩阵 V 由 14 个 向 量 空间 组 成 ， 其 中 VV (A, ) 是 一 个 
216 x 12 的 矩阵， (A) 是 一 个 216 x 6 的 矩阵 ， 
V(B,) 是 一 个 216 x 10 的 矩阵， (B,) 是 一 个 
216 x BAJERE, V (En) -V CES ) 均 为 216 x 18 
的 矩阵。 这 14 个 独立 的 对 称 子 空间 用 于 将 216 x 
216 的 刚度 矩阵 和 质量 矩阵 化 衔 成 14 个 维 数 更 小 
的 子 和 矩阵 。 人 例如， 结构 的 切线 刚度 矩阵 被 矩阵 
V (4，) 转换 成 . 

K (A,) =-V(A,) "+ K-V(A,) = 


23212310 0 0 17535 132867601 0 
0 93630 9961 -25898 0 0 
0 996] 54783 -16186 0 0 


17535 -25898 -16186 21651565 12388557 0 
12388557 49488403 34281 
34281 37037145 


1328673 0 0 
0 0 0 0 


图 7 给 出 了 由 对 称 法 求 得 的 前 100 阶 频率 ， 
并 与 常规 数值 法 计算 结果 进行 了 对 比 。 男 外 ， 为 
了 研究 初始 预 应 力 对 结构 自 振 频率 的 影响 ,我们 
利用 群 论 方法 对 不 同 预 应 力 水 平 的 Ci, HRB 
项 结构 进行 了 参数 分 析 ， 其 中 初始 预 应 力 水 平分 
别 取 为 0、0.25t、0.5t、t、2t 和 4t。 
40 
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图 7 ”不同 预 应 力 水 平 下 Ci PRS HT 100 阶 频率 
由 图 可 知 ， 前 28 阶 自 振 频 率 值 非常 小 ， 且 密 
集 式 分 布 。 对 称 法 所 得 结果 与 常规 数值 法 十 分 吻 





HI H2 N1 N2 C1 C2 C3 
24 36.7 24.9 14.2 66.3 32.8 215 
1637 1307 469 750 —305 —126 -52 


合 ， 其 最 大 误差 为 le-9。 此 外 ， 算 例 再 次 证 明了 
对 称 法 比 常规 数值 法 更 加 高 效 ， 它 们 求 取 广义 特 
征 值 所 消耗 的 时 间 分 别 为 0.063s 和 0.3594s。 与 一 
阶 无 穷 小 机 构 位 移 模 态 相 关联 的 前 14 阶 模 态 的 频 
率 值 大 小 直接 取决 于 初始 预 应 力 的 大 小 ， 如 图 7 
AR. BRAK (2) 可 以 发 现 ， 低 阶 频率 值 与 
预 应 力 大 小 的 平方 根 成 正比 。 事 实 上 ， 初 始 预 应 
力也 较 好 地 改善 了 该 索 杆 结构 的 刚度 ， 提 高 初始 
预 应 力 的 水 平 后 ， 结 构 的 前 50 阶 自 振 频率 均 显 著 
增加 。 

图 8 描绘 了 结构 的 前 12 阶 自 振 振 型 ， 可 以 注 
意 到 ， 索 穿 顶 结构 存在 许多 重复 的 特征 值 和 等 效 
的 特征 向 量 ， 例 如 模 态 1-2、 模 态 3-4、 模 态 $- 
6、 模 态 8-9 以 及 模 态 11-12 。 其 实 ， 这 种 重 根 现 
象 在 对 称 结构 中 普遍 存在 ， 这 是 因为 对 于 1<i< 
pP， 对 称 子 空间 V (En) FIV (E) 的 特 值 征 计算 
是 完全 相同 的 。 另 外 ， 对 称 子 空间 V(E ) 的 振 型 
并 不 能 保持 其 自身 全 对 称 性 (Ci,, ) ， 并 且 可 能 
降 阶 为 更 低 的 对 称 性 ， 壁 如 仅 拥 有 单一 镜像 对 称 
EVE ( 对 称 群 C, ) ， 见 图 8。 图 中 第 7 阶 模 态 的 振 
动 形状 和 运动 趋势 可 以 从 子 空间 4, 得 到 : 根据 对 
称 群 的 降 阶 属性 ， 该 振 型 具有 C ,旋转 对 称 性 。 
另外 ， 第 10 阶 振 型 之 所 以 保持 C ,全 对 称 性 称 ， 
是 因为 它 来 自 于 对 称 子 空间 4, 中。 

4 结论 
本 文 提出 了 一 种 基于 群 论 的 对 称 方法 ， 能 够 
高 效 求解 对 称 型 预 应 力 结构 的 广义 特征 值 问 题 
(如 自 振 分 析 ) 。 根 据 所 推导 的 对 称 坐 标 系 ， 巧 
妙 地 将 结构 的 切线 刚度 矩阵 和 质量 矩阵 进行 了 分 
块 对 角 化 ， 从 而 明显 简化 了 广义 特征 值 问题 的 求 


解 复杂 性 ， 显 著 降 低 了 计算 工作 量 。 针 对 三 种 不 


同 对 称 类 型 的 预 应 力 结构 ( C, 对 称 案 杆 结构 、 
C, ,对称 索 网 结构 、C,, 对 称 Levy 型 索 穹 顶 ) 的 
算 例 均 表明 ， 文 中 所 提出 的 对 称 法 是 准确 有 效 、 
可 行 的 。 
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图 8 C12v 对 称 Levy 型 索 写 顶 结构 的 前 12 阶 自 振 振 型 
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